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ABSTRACT 

Using observational data, we build numerical N-body, hydrodynamical and com- 
bined equilibrium models for the spiral galaxy NGC 5247. The models turn out to be 
unstable towards spiral structure formation. We simulate scenarios of spiral structure 
formation for different sets of equilibrium rotation curves, radial velocity dispersion 
profiles and disk thickness and demonstrate that in all cases a simulated spiral pattern 
qualitatively agrees with the observed morphology of NGC 5247. We also demonstrate 
that an admixture of a gaseous component with mass of about a few percent of the 
total mass of the disk increases a lifetime of a spiral pattern by approximately 30%. 
The simulated spiral pattern in this case lasts for about 3 Gyr from the beginning of 
the growth of perturbations. 
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1 INTRODUCTION 

Attempts to understand the phenomenon of spiral 
structure in galaxies have a long history. It has become 
clear nowadays that the spiral structure is a density wave 
(or waves) propagating in a multi-component stellar-gaseous 
disk. A universal mechanism for the generation of such spi- 
ral density waves that successfully explains the rich morpho- 
logical variety of spiral galaxies does not yet exist. To de- 
scribe the observed patterns, a few spiral- generation mech- 
anisms are usually invoked. Some researchers treat the spi- 
ral structure as long-lived global modes tha t last in the 
galactic disks for tens o f galactic rotations ("B ertin et. al.l 
Il989al lbl: iBertin and Lh] [1996^. Others consider the spiral 
structure as a transient phenomenon, so that the spiral 
pattern c hanges many t imes during the galactic evolution 

i see, e .g. JSellwoodll201lL for references) . iGerola and SeidenI 
19781 ) even suggested that spiral structure is merely a prod- 
uct of recent star formation, outlined by the new-born stars. 
On observational ground, the probl em of the origin of sp i- 
ral structure was raised recently by Eskridge et. al.] (|2002[ ): 
lElmegreen et all (|201lh : lKendall et aD (l201lh . These authors 
found that galaxies which are optically grand design spirals. 



show a grand design structure in the near-infrared as well. 
The spiral structure therefore is not merely a manifesta- 
tion of a recent star formation but represents also a spiral 
structure in the underlying ste llar mass distribution. The 
paper bv lElmegreen et "ail dSml i shows that some flocculent 
spirals remain flocculent even in the IR. This observational 
evidence rules out at least one approach which suggested 
that the spiral structure in galaxies is a pattern of a recent 
star formation and does not involve an underlying stellar 
mass distribution. To understand to what extent the vari- 
ous suggested explanations are consistent with the observed 
phenomena of spiral structures, a detailed, quantitative com- 
parison between the theoretical predictions and the observed 
properties of the spiral patterns in galaxies should be under- 
taken. In this paper, we continue our efforts to compare the 
properties of the observed spiral structure with the theoret- 
ical p redictions. In previous papers (Korchagi n et al.ll2000l , 
120051 ). we made such a comparison using a hydrodynami- 
cal approach to model the dynamics of galactic disks. In 
this paper, we simulate three-dimensional collisionless, and 
collisionless-gaseous models of galactic disks, aiming to de- 
termine the model parameters that most closely match the 
observed spiral pattern in the galaxy NGC 5247. We choose 



2 Khoperskov et al. 




Figure 1. Lef t : opt ical 
(lEskridge et. al.l (I2OO2I )) 



image of the NGC 5247 galaxy 
Right: HII regions outlining the spiral 



this galaxy for our comparison because it is well studied ob- 
servationally, and its rotation curve as well as disk velocity 
dispersion and luminosity distribution have been accurately 
determined. 

In the cases that have been studied, the admixture of a 
:aseous c omponent increases the lifetime of a spiral pattern 
Semelin and Combesll2QQ2l : lFuxlll999l V In general however, 
the lifetime of the spiral patterns in isolated galaxies and 
the role of gas remain an open issue, that we address in this 
paper. 

In Section 2 we discuss the observational properties of 
the galactic disk of NGC 5247. In Section 3 we present the 
basic equations and discuss the disk equilibrium. Section 4 
presents the results of numerical simulations of a gaseous 
and a collisionless stellar model of the galaxy NGC 5247. 
Section 5 presents the results of our simulations of stellar- 
gaseous models, and in section 6 we present the concluding 
remarks. 
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Figure 2. Radial dependence of the observed velocity disper- 
sion of the stellar disk of NGC 5247 . Filled circles — data 
from Ivan der Kruit a nd Freernanl (|l986l ) , open squares — mea- 
surements of lsottemal (|l993l ). 



I Zhao et aD (|2006l l tackled this problem by using the Jeans 
equations applied along the vertical axis of the disk, and 
estimated the disk thickness of NGC 5247 to be zo — 1.5=b0.6 
kpc. With this value, the ratio of the radial to vertical scale 
length is r^/zo = 3.2tQ9. 

The line-of-sight velocity dispersion as a func- 
tion of galactocentric radius for NGC 5247 ha s also 
been determined from observations (|Bottemal Il993l : 
Ivan der Kruit and Freeman| [l986, see Fig. (2). The measure- 
ments by both groups agree with each other and show ap- 
proximately exponential decrease of the velocity dispersion 
with radius, which is typical for disk galaxies. 



2 NGC 5247: CONSTRAINTS FROM 
OBSERVATIONAL DATA 

NGC 5247 shown in Figure [1] is a nearby spiral galaxy 
which probably belongs to the Virgo Supercluster of galax- 
ies (Bingelli et a. 1985). It has a strong, well devel- 
oped two-armed "grand-design" spiral pattern and does 
not demonstrate any signs of interactions with other galax- 
ies (|Considere and Athanassoula 1988 ). The spiral pattern 
however looks somewhat lopsided. A low-amplitude spiral 
pattern is detected in an interarm region. However the mor- 
phology of the galaxy is dominated by the two-armed spiral. 
The photometric data do not point to any noticeable bar in 
the central regions of the disk, therefore NGC 5247 is a clas- 
sical example of a grand-design spiral galaxy. 

Our numerical raodelin g is based on photometric data 
taken by I Zhao et aP (l2006l V We take as the radial scale 
length of the disk Vd — 4.8 kpc as determined by these au- 
thors. The inclination angle of NGC 5 247 is rather unc er- 
tain. Me asurements range from i = 20° (IZhao et al.ll2006l ) to 
i — 40° (HyperLeda") . In this paper we adopt an inclination 
of 28° which is close to an average from these references, 
and a distance to the galaxy of 17.4 Mpc (|HvperLedah 

The vertical scale of the disk (zo) significantly influences 
its stability properties. Unfortunately, an accurate estimate 
of the vertical thickness of the galactic disk is not simple. 



3 MODELING NGC 5247: BASIC EQUATIONS 

We assume that the potentials of the dark halo and 
the stellar bulge are static and remain axisymmetric. Such 
an assumption does not affect the dynamics of a galactic 
disk because a contribution of high velocity dispersion com- 
ponents to the spiral density waves is insignificant even if 
they involve most of the galactic mass (Marochnik 19^). 
The galaxy does not have a large central bar, and we ex- 
clude from a consideration the models with a bar which is 
conside red to be an effective generator of the spiral densit y 
waves ("Korchagin and Shevelev" ('l98l'); B uta et all (|2009l l: 
K aufmann and Contopoulos (1996)). Second, observational 
data do not allow us to estimate the degree of asymmetry 
of spheroidal subsystems, and the emergence of new free 
parameters greatly complicates the analysis. Moreover, an 
interaction with a live bulge or a halo can not be decisive. 

To simulate the dynamics, we use a set of equations 
which self-consistently describe the behavior of a stellar- 
gaseous disk in equilibrium with an external gravitational 
field ^ext- The set of equations for the stellar disk consist- 
ing of N gravitationally interacting particles is written as: 

^ = -V[*.(r,) + *e.t(ri) + *p(r,)] , i = \,...,N. (1) 

Here is the gravitational potential of the stellar compo- 
nent of the disk, is the potential caused by the gaseous 
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component, and ^ecct is the combined gravitational poten- 
tial of the galactic halo (^^) and the bulge (^b)- Therefore 

The set of equations describing the dynamics of the 
gaseous component are given by: 



volume approximation of the variables is used: 



+ V-(p,u) = 0, (2) 

+ U- V(PpU) + Vp + PgV [^g + ^ext + ^s] = , (3) 



dt 

dt 
OE 

— + V- [{E + P)U] + PgUV [^g + ^ext + ^s] = , (4) 

here pg, p, and u are the volume density, the pressure and 
the velocity vector of gas. The energy density of gas E is 
given by the expression E = p^|u|^/2 +p/(7 — 1). We use 
an adiabatic equation of state pg = KsP^ for the gas, where 
the value of the adiabatic index is assumed to be 7 = 1.1 
and the constant Ks was chosen so that the sound speed is 
equal to Cs = 10 km s~^ in central regions of the disk. The 
sound speed of the gas decreases slowly with radius, reaching 
approximately 8 km s~^ at the periphery of the disk. Self- 
gravity is taken into account by the Poisson equation applied 
for the gaseous and stellar components: 



(5) 
(6) 



To model the gravity of the halo, we use the potential: 

(7) 



1, 



Here ^ = r/a^, is the scale length of the halo potential, 
Mh is the mass of the halo within radius rh — 20 kpc and 
Ch = GLh{Th/oLh — arctg(rh,/ah,) ) • Such a choice provides a 
constant rotation velocity at large radii r > 2ah^ so we have 
a rotation curve typical for disk galaxies. 

The King model with a cutoff a density distribution of 
large radii r > is used to model the potential of a stellar 
bulge: 



^6 = 



In 



+ A 1 + 



(8) 



here Vh is the bulge scale length, Mh is the bulge mass within 
rr and ^ = In ( + Jl+f^) ) - ^ ^"Z"" . 

The observed rotation curve of the disk can be repro- 
duced by varying the parameters of the dark halo (M^,, a^), 
the stellar bulge (Mb, Vh) and the stellar disk (Ms, r^, zq). 

To analyze a behavior of a one-component collisionless 
model, we assume that the gas density pg is equal t o zero, 
and solve equation (|6} using TREE-code (Barnes an d HutI 
Il986l l adapted to the parallel calculations. The number of 
particles used in our collisionless models range from A/" = 10^ 
to 10^. To build a self-consistent model of a stellar-gaseous 
disk, the set of hydrodynamical equations for a gaseous disk 
rotating in a self- consistent gravitational potential is added 
to the equations of motion for a stellar disk. 

We apply a TV D MUSCL type scheme in the cylindrical 
coordinate system (|van Leerl [l979) to solve the equations 
([2])-(|3|). To preserve conservation laws at a grid level, a finite- 



i+1/2 



: q^ + 0.25((/c + 1)L>+ + (/c - l)D-))^ , 



here {q) 1X0)1+1/2 are the conserved variables at a bor- 
der of i-ih cell and at its center correspondingly; D+ — 
minmod{d+,bd-), D- — minmod{d-,bd+), where d+ 
are the variations of the conserved variables to the left and to 
the right of a cell border. The parameter 6 = (3 — k)/{l — k) 
determines an order of approximation in space. We use the 
value k = 1/3, corresponding to the third order of spatial 
approximation. The function 

, sign(x) + sign(y) /i , , ix 
minmod[x,y) = rmn{\x\, \y\) 

is used to reconstruct the discontinuous functions between 
the nodes of the computational grids. 

Interaction between the stellar and the gaseous com- 
ponents of the disk occurs due to gravity. Similar to the 
one-component models, such interactions are computed us- 
ing the TREE-code. 

At the beginning of the simulations, the stellar disk is 
in an equilibrium in the radial and the vertical directions 
with its density given by: 

Ps = Pso exp(-r/rd) • A{z/zo) , 

here r^, zq are radial and vertical scales of the stellar disk. 
Function A{z) determines a vertical distribution of density 
of the stellar disk: A(z) — c]i~'^ {z / zq) . 

The equilibrium of the stellar disk in the vertical direc- 
tion is determined by solving the Jeans equation: 

d f 2dps\ 2 [ dps\'^ , 



^2 2 
2 d 



7^ + 4^^/^' {Ps+E + pext (^)) + ^ J = ' 



dz'^ 



d{rpsarz 



AnGr dr r dr 

Here Cz is the stellar velocity dispersion in the vertical di- 
rection, arz = {uw) is the result of averaging the product of 
the radial u and vertical w velocity components, pext is the 
total density of the halo and the bulge. 

By using Jeans equation in the radial direction, we de- 
termine the rotational velocity of the stellar disk in the disc 
plane z = (jKhoperskov et aLllioToh 



K' + c?{l- J 



r d{psCr) r darz 



psCr 



dr 



Cr dz 



(10) 



here Vc is the circular velocity of a test particle in an ax- 
isymmetric field ^: /r = , Cr/c^ = 2^/^, Q = V/r 

and — Y^^^^Ml + ^^^V(2^^^0) i^ ^^e epicyclic frequency. 

The gaseous disk extends in our models beyond the op- 
tical radius of the stellar disk to the distance r = IQra- The 
density distribution of the gaseous disk is modeled by the 
expression 

P9 = pAI - r/(16rd)) • B{z) 

where the function B[z) is determined by the balance of 
the gravity force and the gas pressure gradient. With such a 
distribution, the gas density decreases by a few times within 
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the optical radius, and goes to zero in the ghost zones of 
the computational grid at r = l^Vd. The Euler equation 
uniquely determines the rotational speed of the gaseous disk: 



300 



— = / - ^ + + 

r pg or or 



(11) 



We build a set of models for a stellar-, and a stellar- 
gaseous disk of the galaxy NGC 5247 such that the model 
predicted quantities like rotation curve and velocity disper- 
sion are consistent with the observational data (and the un- 
certainties associated with the observations). In the follow- 
ing sections, we explore the model parameter space and dis- 
cuss the dependence of the morphology of spiral pattern on 
the rotation curve, the stellar velocity dispersion, the verti- 
cal thickness of a stellar disk, and mass ratio of gaseous to 
stellar component of the galactic disk. 



4 DISK EQUILIBRIUM 

4.1 Rotation curve 

A small inclination angle of the galaxy influences sig- 
nificantly a reconstru ction of its r o tation curve. Using the 
Tully-Fisher relation, IPatsis et al.l (|l997f ) found the maxi- 
mum rotational velocity of the galaxy to be: Vmax = 205 
km/s. Other measurements give the values Vmax = 165 ± 20 
km/s (IBotte ma 1993) and Vm ax — 3 00 km/s. Following the 
paper of iContopoulos and Grosbo 3 ([19861 we build a set 
of rotation curves that have a quasi-solid rotation in the 
central regions of the disk, and are fiat in its outer regions. 
Such a shape of the rotation curve is typical for a large spiral 
galaxy. To build more realistic models, we take into account 
the gravitational potentials of the halo, of the bulge and 
of the galactic disk. The parameters of the potentials and 
corresponding properties of the rotation curves are listed in 
Table [11 To take into account observational uncertainties in 
our knowledge of the rotational curve of the galaxy, we vary 
the rotational velocity of the disk. Namely, the model A has 
the largest rotational velocity, B - middle one and model 
C has the lowest rotational velocity. The resulting rotation 
curves are shown in Fig. (3] A large spread in the rotational 
velocities reflects an uncertainty of observational data for 
this galaxy. 



4.2 Velocity dispersion 

To constrain the radial component of the velocity dis- 
persion of the stellar disk Cr we use observed line-of-sight 
velocity dispersion data of NGC 5247. For the nearly face- 
on orientation of NGC 5247, the line-of-sight velocity dis- 
persion is close to the component of the velocity dispersion 
in the direction perpendicular to the disk. 

To mimic the observational uncertainties in our knowl- 
edge of the velocity dispersion, we choose three profiles as 
shown in Fig. [U In this Figure, index 1 corresponds to the 
hottest disk, and index 3 corresponds to a model with the 
lowest velocity dispersion. 

Observational data of the external galaxies show that 
the ratio of the vertical velocity dispersion to the radial ve- 
locity dispersion is approximately constant along the radius 
of a galactic disk, and varies from galaxy to galaxy within 
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Figure 3. Rotation curves for one-component models. The maxi- 
mum values of the rotation curves are equal to 245 km/sec (Model 
A), 205 km/sec (Model B) and 153 km/sec (Model C). 
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Figure 4. Radial velocity dispersion as a function of galactocen- 
tric distance. In model 1 — Cz/cr = 0.43, model 2 — Cz/cr = 0.6, 
model 3 — c^/cr = 0.8 



Czlcr = 0.3 - 0.8 (iKorchagin et al.ll2000l : iKhoperskov et al.1 
I2OIOI ). In our models we choose this ratio equal to Cz/cr — 
0.43; 0.6; 0.8 in central regions of galaxy. 



4.3 Disk thickness 

A direct measurement o f the vertic a l scale of the disk 
of NGC 5247 is impossible. IZhao et al.1 (|2006l ) solved this 
problem by applying the Jeans equation in the direction 
perpendicular to the collisionless stellar disk. They found 
the value of the disk vertical scale to be zq — 1.5 ± 0.6 
kpc. Taking into account this result, we additionally consider 
two variations of the model C3 with larger thickness of the 
disk equal to zq — 2.1 kpc, and the lower disk scale length 
ZQ = 0.9 (model D3). 

A variation of the thickness of the stellar disk is taken 
into account in a set of models D with thickness of the disk 
decreasing from model Dl to model D3. 
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5 ONE-COMPONENT MODELS 



5.1 Hydrodynamical model 



The simplest model to study the dynamics of an un- 
stable galactic disk is the hydrodynamical approximation. 
In this approximation, a stellar-gaseous gravitating disk is 
represented by a gravitating compressible fluid. There are 
some arguments justifying the hydrodynamic approx imation 
for a descri ption of stellar disks: .Marochnik (J^^) , iHunterl 
(|l979l ). and ISvgnet et al.1 (|l987l l demonstrate that the be- 
havior of perturbations in collisionless disks can be approx- 
imately described by introducing an isotropic pressure with 
polytro pic constant 7 = 2. Us in^ the hydrodynamic approx- 
imation iKikuchi et aP (|l997l ) made a direct comparison of 
the global stability properties of gravitating disks with the 
solution of t he collisionless Boltzman n equat ion given by 
IVauterin and Deionghe (1996). Kikuc hi et al.1 ^997) found 
a qualitative, and in most cases a good quantitative agree- 
ment between the results found in a hydrodynamic approx- 
imation and by the direct solution of the Boltzmann equa- 
tion. 

To study the dynamics of the disk of NGC 5247, we use 
the parameters of model CI, listed in Table [Tl This model 
was studied by using the linear global stability analysis as 
well as two-dimensional nonlinear simulatio ns. Both meth- 
ods are des cribed in detail i n pape rs by Korchagi n et al.l 
([2OOQ ) and iKorchagin eraP (|2005l l. The linear analysis 
yields the most unstable global mode as the m = 2 mode 
with eigenvalues Iie(uj2) — 1.67, lm.{ijj2) — 2.77. Its near- 
est competitor, m = 3 mode, has the imaginary part of 
the eigenfrequency that determines the mode's growth rate 
equal to lm.{uj2,) — 2.43. With such values, the dynamics of 
perturbations is determined by a two-armed spiral pattern. 

To analyze quantitatively the large-scale morphology 
of the spiral structure formed in the disk, we use a 
Fourier decomposition of the surface density S(r, = 
Qs{r, z)dz of the disk. 

Figure [5] plots the time dependence of the global Fourier 
amplitudes for modes m = 1 — 4 



ir lo ^ir,<p)rdrd<p 



(12) 



growing from random initial density perturbations. Snap- 
shots of density perturbations in the nonlinear simulations 
are shown in Figure {|6)), which clearly illustrates how the 
two-armed spiral structure evolves out of the initial random 
density perturbation. 



5.2 Collisionless models 

Although the hydrodynamical approximation repro- 
duces qualitatively the morphology of NGC 5247, the life- 
time of the spiral pattern is short. After approximately 0.3 
Gyrs the spiral pattern is destroyed in these hydrodynam- 
ical models and they no longer resemble the observed spi- 
ral structure. A mechanism that destroys the spiral struc- 
ture is not completely understood. Using hydrodynamical 
model, iLaughhn etld] (|l997l l showed that a nonlinear spi- 
ral mode causes a hollowing out of the surface density profile 
in a vicinity of a corotation resonance. This stops the spiral 
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Figure 5. Time evolution of the Xo^^^i^Am) for m = 2 — 4 global 
Fourier amplitudes. The most unstable mode is the m = 2 spiral. 
The instability saturates at level of log;|^Q(Am) = 0.5. 
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Figure 6. Time evolution of density perturbations in galactic 
disk in model CI. The spiral is in a counterclockwise rotation 



growth and eventually destroys the spiral itself. For the col- 
lisionless disks, the mechanism remains unclear. We explore 
therefore the coUisionless models to study to what extent the 
collisionless nature of the galactic disk increases the lifetime 
of the spirals. 

Using distributions shown in Figure [7] we build a set 
of models for a stellar disk that have stability Q-parameter 
larger then unity: 



200 



1000 




8 12 
distance, [kpc] 

Figure 7. Initial distributions of the stellar-gaseous disk param- 
eters in model E2: velocity of stars (open circles), velocity of gas 
(fill circles), radial velocity dispersion of stellar disk (fill trian- 
gles), speed of sound in gaseous component (open triangles), stel- 
lar surface density (fill squares) and gaseous surface density (open 
squares) 
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Figure 9. Time dependence of angular velocity of two- armed 
pattern in model C3. 
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(13) 



If condition (|13p does not hold for the whole disk, then 
the perturbations grow significantly during a short evolution 
time, and the parameters of the disk change within 1-2 disk 
rotations. We exclude such cases from our consideration. 

The rotation curve significantly influences the dynam- 
ics of a collisionless disk. Figure [8] shows evolution of the 
Fourier harmonics Am in three models with different rota- 
tion curves. In the fast rotating model A3, the disk is stable 
due to the presence of a massive halo. The model with the 
fastest disk rotation therefore can not explain the observed 
spiral pattern in the galaxy NGC 5247. In model B3 (mid- 
dle frame) , the perturbations slowly grow with the most un- 
stable mode two-armed spiral. During approximately three 
Gyr from the beginning of simulations, the wave amplitude 
reaches about 2-5% of disk unperturbed values which is sig- 
nificantly lower then the amplitude of the observed spiral 
pattern. Apparently, the models A and B can not explain 
the observed spiral structure of the galaxy NGC 5247. 

Let us examine model C, which has a more massive 
stellar disk compared to models A and B. The minimum 
value of the Toomre Q-parameter indicates that this model 
is gravitationally unstable. Similar to model B, there is a 
well defined two-arm global spiral structure that extends to 
about four scalelengths in the disk. The amplitude of the 
perturbations is much higher compared to the other cases. 
During the initial stages (t < 0.5 x 10^ years) one can see 
a three-arm spiral, and the spiral pattern is a superposition 
of a two- and three-armed modes. However, in the nonlinear 
stage a two-armed pattern dominates. Figure [9] shows a typ- 
ical dependence of an angular velocity of a two-armed spi- 
ral pattern generated in model C3. Between approximately 
0.2 — 0.7 Gyr when an exponential growth of a spiral per- 
turbation is observed, a spiral perturbation has an angu- 
lar velocity equal approximately to 22 km s~^ kpc~^. After 
0.8 Gyr, a nonlinear saturation starts, the angular velocity 
decreases to 18 km s~^ kpc~^ and remains approximately 
constant during the rest of simulations. 

Despite significant difference in the values of radial ve- 
locity dispersion used in our numerical simulations, the spi- 
ral patterns are quite similar as can be seen by compari- 
son models CI, C2, and C3 in Figure [TTl On a qualitative 



t = 0.5 Gyr 



t = 1 Gyr 



t= 1.5 Gyr 




Figure 10. Evolution of the surface density in the numerical 
simulations of the collisionless disks with different rotation curves 



level one can conclude that morphological properties of spi- 
ral structure are independent on variation of the radial veloc- 
ity dispersion within observational errors. The growth rates 
however are different ranging from 0.7 Gyr in the coldest 
disk to 1.5 Gyr for model with the largest velocity disper- 
sion (Figure [12]). 

The value of the disk vertical scale used in our simu- 
lations is close to the value estimated by Zhao et alj (|2004[ ) 
1.1 — 2.5 kpc, implying a large range of values to be explored. 
Figure [13] shows that the spiral patterns in collisionless cases 
D1-D3 are similar, and independent of the disk thickness. 

Taking into account the third dimension, i.e. the disk 
thickness, is important in disk dynamics to simulate cor- 
rectly the density growth in spiral arms. The vertical mo- 
tions lead to a decrease of the growth rate of an unstable 
spiral wave. Under similar conditions, the amplitude of the 
spiral wave grows faster in a thinner disk (Fig. [13]). In the 
model with a relatively large scale length of zq = 2.1 kpc, 
the structure of the spiral arms and their lifetimes are qual- 
itatively similar to the behavior of model C3. However, the 
disk with a small thickness is significantly more unstable, re- 
sulting in a rapid growth of multi-armed perturbations and 
formation of a flocculent spiral pattern. 

In general, the lifetime of the global two-armed spiral 
pattern in our collisionless models is about 1.5 Gyr. As was 
mentioned, a gaseous component in the disk increases the 
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lifetime of the spiral structure. In the following section we 
will consider this effect in detail. 



6 TWO-COMPONENT MODELS 

In our two- component models, we take into account 
the gravitational interaction between the gaseous and the 
stellar components of the disk. Effects of star- format ion or 
exchange of matter between the stellar and gaseous sub- 
systems are neglected. It is known that gas has a strong 
destabilizing effect on a gravitating disk even if its surface 
density is much smaller than the surface density of the stel- 
lar disk. We build our stellar-gaseous models based on the 
one-component model C3, and by varying the total gas mass 
within the radius of the stellar disk. Namely, we analyze the 
models with the relative mass of a gaseous component equal 
to: 

Mg/Ms = {0.01; 0.05; 0.1; 0.2} . 

For a two-component gravitating disk, each subsystem 
has its own spatial density distribution, rotation curve and 
velocity dispersion profile. Figure [71 shows initial equilibrium 
distributions of basic quantities for a stellar-gaseous disk. A 
difference between the rotation curves of stellar and gaseous 
components is caused by a difference in their radial veloc- 
ity dispersions, and influence of the velocity dispersions on 
disk equilibrium. The stellar disk has an exponential den- 
sity distribution with a radial scale length of = 4.8 kpc. 
The velocity dispersion of the stellar disk in the direction 



perpendicular to the disk decreases from 70 km s~^ in its 
central regions to about 20 km s~^ at disk periphery. The 
gaseous disk is cold compared to the stellar one. Its velocity 
dispersion is close to 10 km s~^ in the disk's central regions 
and changes slowly within the stellar disk. 

For a two component gravitating disk, the stability cri- 
terion needs to be modified. A few such modifications of 
the stability parameter have been suggested in the liter- 
ature. A simple stab ility criterion has been proposed by 
IWang and SiS (Il994l l: 

111 

Here Qs and Qg are the values of the Toomre Q-parameter 
for the stellar and t he ga seous components respectively. 
"Rome o and WiegertI (|201lh demonstrated that such an 
approach is not su i table for three-dimensional systems. 
iRomeo and WiegertI (|201l[ ) suggested an expression for the 
effective stability criterion that takes into account the finite 
thicknesses of a stellar and a gaseous disk: 

{W 1 
rri ^ ^ rri ^ 5 TgQ S ^ TgQ g] 
+ TO >TO 

Here W = 2c,c,/(c^ + c?), = 0.7+0.8{c, /cr)s, Tg = 0.7+ 
{cz/cr)g. The radial dependence of the stability parameters 
Qs, Qg and Qrw are shown in Figure [Ml 

The evolution of the surface density perturbations in 
a stellar disk in models with different admixtures of gas 
is shown in Fig. [151 During t < 1 Gyr, a spiral pattern 
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Figure 11. Evolution of the perturbed surface density of the 
stellar disk in simulations with different radial velocity dispersion 
profiles. 

is formed in all models. Later on (t > 1.5 Gyr), the spi- 
ral wave dissipates in the purely collisionless model, and 
becomes indistinguishable from a background of low-scale 
density perturbations. A small gas admixture with a mass 
of a few percent of the total mass of the disk significantly 
increases the life time of the spiral structure, and spiral pat- 
tern exists up to 3 Gyr. If mass of the gas is increased to 
one-tenth of the total mass of the disk, then a strong gravi- 
tational instability is created, leading to the formation of a 
complex, highly transient multi-armed spiral structure. The 
disk changes significantly from its initial state during about 
0.5 Gyr. In such a case, star- format ion should be taken into 
account which is beyond the scope of our paper. 

The self- consistent stellar-gaseous models allow repro- 
ducing a long- livin g two- armed spiral patter n . We find 
however, similar to [Elmegreen and Thomasson (logs') and 
iThomasson et al.] ([l990), that the morphology of the spiral 
pattern evolves with time, with overall lifetime of the spi- 
rals being about ten disk rotations. In comparison, lifetime 
of spirals in the one-component models do not exceed 1.5 
Gyr. 

Figure [16] shows the azimuthal density distribution at 
r = 4.9 kpc in the gaseous, and in the stellar components. 
The figure clearly demonstrates formation of shock fronts 
in the gaseous component located at the inner edge of the 
stellar arms. 



t = 0.5 Gyr t = 1 Gyr t= 1.5 Gyr 




:..:?.^-f^^: R , . , : , . . . , nrirr^'?^. 
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Figure 13. Evolution of the surface density perturbations in stel- 
lar disk in numerical simulations with different thickness of the 
stellar disk: B3 - zq = 0.9 kpc, B2 - zq = 1.5 kpc, Bl - zq = 2.2 
kpc. 




0.5 r ■ . ■ ■ ■ 1 

4 8 12 16 20 
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Figure 14. Radial distribution of Toomre parameter: gaseous, 
stellar and both stellar-gaseous disk for model E2 
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Figure 15. Comparison of the evolution of the surface density 
perturbations of a stellar disk in numerical simulations with dif- 
ferent amount of a gas-component 

The shock front is formed when the supersonic flow of 
inter stellar gas pass es through the density wave of the spiral 
arm (|Robertsl ll969'). Distribution of gas along the azimuth 
indicates the position of the shock wave on the inner edge of 
the stellar density wave in our simulations. In the presence of 
a strong shear flow the syste m of shock fronts is unstable to 
hydr odynamical instabilities (|Wadall20o3 : iKhoperskov et all 
l201lh . We can not reveal such an instability due to the low 
resolution of our simulations. 

Figure [T71 shows azimuthal dependence of the stellar 
surface density at different radii in the disk. As one can judge 
from the figure, the amplitude of the spiral perturbation 
decreases with radius in a stellar component. The shift in 
phase of the density maxima with radius reflects the spiral 
nature of the perturbations. 



7 DISCUSSION 

7.1 Toomre's parameter 

ISellwoodI (I 201lh reviewed observational as well as theo- 
retical arguments and came to the conclusion that spiral pat- 
terns are short-lived. A qualitative counterargument against 
this picture is that if spirals recur often enough a large frac- 
tion of galaxies should have a low-amplitude spirals. This is 
not the case in the most of observed cases. Below we discuss 



10^ 




azimuth angle 

Figure 16. Azimuthal layer of the gaseous and stellar surface 
density at t = 1.5 Gyr and r = 4.9 kpc for model E2 




azimuth angle 



Figure 17. Azimuthal layers of a stellar-component surface den- 
sity for various values of radius — from top to bottom: r = 
{2.88; 4.8; 6.72; 8.64} kpc for model E2. 

in detail two factors, considered by Sellwood as arguments 
against long-lived spiral structure. 

One of his arguments is that the theory of spiral modes 
requires galactic disks to be dynamically cool. Namely, 
Toomre's axisymmetric stability parameter should range 
within 1^(3^1.2. Asa test of applicability of a modal ap- 
proach in such systems, Sellwood simulated the dynamics of 
a disk that has Q-parameter equal to unity everywhere ex- 
cept the disk's central regions where the Q-parameter rises 
steeply. Via this experiment Sellwood demonstrated that 
such a low value of Q-parameter supports vigorous collec- 
tive responses that change the disk properties during a short 
dynamical time. 

We repeat the experiment of Sellwood in 3D simulations 
assuming that the disk of the galaxy NGC 5247 has Q = 1 
everywhere except its central regions. We find a q ualitative 
agreement of our simulations with the results of ISellwoodI 
During approximately one disk rotation a complex 
multi-armed transient spiral structure is formed that causes 
disk heating and a quick dissipation of the spiral waves. 
Figure [18] shows the radial profile of disk radial velocity dis- 
persion at different moments of time. After a violent disk 
heating, it comes to a stationary stage. 
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Figure 18. Radial distribution of the radial velocity dispersion 
at different times: dotted line is initial statement {t = 0), dash- 
dot line is distribution after one rotation period (t = T), full line 
is stationary distribution (t/T = 3 ^ 10). 

The assumption about a radial profile of Toomre's Q- 
parameter invoked by Sellwood is not related to the spiral 
galaxies. It is known from the direct measurements that ve- 
locity dispersion in the galactic disks decreases exponentially 
with radius and has a larger scale length compared to the 
scale length of the surface brightness (or density) . In combi- 
nation with the fact that the disks are nearly flat and self- 
gravitating in the vertical direction, this leads to a bell-like 
shape of Toomre's Q-parameter rising towards the center of 
a disk and on its periphery. With such a Q-profile, the disks 
are globally unstable and support a spiral structure during 
tens of galactic rotations. 

We note also that if real galaxies were to have a Q- 
parameter profile close to that assumed by Sellwood, they 
would be in a phase of violent evolution and have a very 
short dynamical time. This would be be difficult to recon- 
sile with the number of observed spiral galaxies and their 
lifetimes. The density wave theory (the disk's instability 
towards the global spiral modes) does not require the Q- 
parameter to be in a narrow range 1 ^ Q ^ 1.2 as it was 
demonstrated by a number of ex amples in hydrodynamical 
simulations by korchagin et al.l (I2O OO. 2005) and in colli- 
sionless models bv iKhoperskov et al.l (120071) . 

7.2 Resonances 

After the work of lToomr3 (|l98lh , the swing amplifica- 
tion of perturbations in a differentially rotating gravitating 
disk is regarded as one of the key mechanisms that amplifies 
the spirals. The swing amplification works in the presence of 
a shearing disk. The orbital clock is faster at the interiors of 
the galactic disks, and the mechanism runs more rapidly at 
smaller radii. The important part of the swing amplification 
mechanism is a propagation of the waves in a radial direction 
with a group velocity. A propagating wave reflects from the 
center and/or from a corotation region, changing from trail- 
ing to leading, and simultaneously changing a sign of the 
group velocity and serves, according to Toomre, as "fresh 
grist to a swing amplified mill" to operate continuously in a 
disk. 

In this mechanism, the inner Lindblad resonance is con- 



sidered as the most important in generation of the spiral 
density waves. This resonance must be shielded to prevent 
a fierce damping of the waves. To check the importance of 
the inner Lindblad resonance in generation of spirals, we 
performed an artificial experiment with our hydrodynami- 
cal model imposing an absorption region in the inner regions 
of the disk where all perturbations are zeroed out. This cuts 
the feedback loop and destroys amplification of spiral per- 
turbations. We performed a series of 2D simulations of the 
unstable 2-D gaseous disk with such an artificially imposed 
"absorption ring" where all perturbations are forced to be 
zero and macroscopic characteristics of the disk are set to be 
unperturbed in a smooth way. Fig. [19] shows angular veloc- 
ity, epicyclic frequency and positions of resonances for such 
a disk. As one can see, the model does not harbor the in- 
ner Lindblad resonance, so an absorption ring was imposed 
artificially to mimic its influence. 

Fig. 1201 shows a result of such simulations. On the right 
panel of Fig. [20] is shown a two- armed spiral pattern devel- 
oping in an unstable 'template" gravitating disk. The left 
panel shows the results of simulations when an absorption 
ring has been imposed close to the disk center. The disk 
remains unstable and a two-armed spiral is developing out 
of noise perturbations. We find that the spiral perturba- 
tions develop in the disk independently of the position of 
the absorption ring, be it inside, or outside the corotation 
resonance. 

An instructive experiment is to impose the absorption 
ring at corotation. The behavior of perturbations changes 
dramatically in this case (Fig. [2T]) . The growth of a two- 
armed spiral is totally suppressed, and a multi- armed spiral 
starts to grow with a four-armed spiral prevailing. These 
experiments illustrate that the qualitative picture of spiral 
growth by swing amplification is not necessarily as relevant 
as once thought. Instead, the corotation resonance may be 
more important in spiral generation. 

The issue of transient spirals is c urrently a topic of 
heated debate. IZhand (|l996l . Il998l . Il999l ) finds that galaxies 
undergo a secular evolution that leads to a redistribution of 
the disk matter and heats the disk stars so they gradually 
rize above the galactic plane in the disk central regions. Both 
effects are the natural ingredients of a presence of the spiral 
density wave. Stars migrate inside and outside the corotation 
region, and a more and more centrally concentrated density 
distribution is achieved with time together with the buildup 
of of an extended outer envelope. Similar effect was fou nd 
in hydrodynamical simulations bv lLaughlin et aT 



8 CONCLUSION 

We have built one-component collisionless and hydrody- 
namical models as well as multi-component stellar-gaseous 
models of the disk galaxy NGC 5247. The models are un- 
stable to m=2 global modes that form a two-armed spiral 
pattern with amplitudes of about 10-20%, qualitatively sim- 
ilar to the observations. 

Our simulations show that in purely collisionless mod- 
els the lifetime of the spiral structure is about a few galactic 
rotations, and does not last longer than one Gyr. An admix- 
ture of a gaseous component with mass of a few percent of 
the mass of the stellar component significantly increases the 
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Figure 19. Angular velocity, epicyclic frequency, and the positions of principal resonances for a gaseous disk (re - corotation radius, 
Voir is the outer Lindblad resonance). The inner Lindblad resonance is not achieved in this model. The units of length and frequency are 
2 kpc and 72 km s"-*^ kpc"-*^ correspondingly. 




Figure 20. The right panel shows a two-armed spiral pattern developing in an unstable "template" gravitating disk (see section 7.1 for 
details). The location of the corotation radius is shown by red circle in both the right and the left panels. The left panel shows the results 
of a simulation when an absorption ring (whose inner and outer extent are shown by yellow circles) has been imposed close to the disk 
center. Note that the development of the 2-armed spiral pattern is not inhibited by the presence of the absorption ring, in contradiction 
to the expectations from standard swing amplification mechanism. 



lifetime of the spiral structure. In the simulated model of 
galaxy NGC 5247, the spiral lifetime increases up to 3 Gyr, 
and comparable with the lifetime of galaxies. 
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circles) has been imposed at the corotation radius. Note that unlike Fig. [20] growth of a two-armed spiral is completely suppressed due 
to the presence of the absorption ring, and a multi-armed spiral starts to form. 
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Table 1. The parameters of one- component (N-body) and two-component models, c^o is the radial velocity 
dispersion in center of the disk, Vmax — the maximum value of the rotation velocity, percentage of gas - 
Mg/Ms, one-component model if equal zero. 



Model name 


Cz 1 Cr 


Cro (km/s) 


ZQ (kpc) 


Vrna. (klll/s) 


Mg/Ms 


A 1 


n /I Q 
U.4o 


1 on 


i.O 




U 


A2 


0.6 


92 


1.5 


245 





A3 


0.8 


70 


1.5 


245 





Bl 


0.43 


120 


1.5 


200 





B2 


0.6 


92 


1.5 


200 





B3 


0.8 


70 


1.5 


200 





CI 


0.43 


120 


1.5 


153 





C2 


0.6 


92 


1.5 


153 





C3 


0.8 


70 


1.5 


153 





Dl 


0.8 


70 


2.2 


153 





D2 


0.8 


70 


1.5 


153 





D3 


0.8 


70 


0.9 


153 





El 


0.8 


70 


1.5 


153 


0.01 


E2 


0.8 


70 


1.5 


153 


0.05 


E3 


0.8 


70 


1.5 


153 


0.1 


E4 


0.8 


70 


1.5 


153 


0.2 



